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Zhang recently conjectured an approximate 50(5) symmetry relating antiferromagnetic and superconducting states 
in high-Tc cuprates. Here, an exact SO{5) symmetry is implemented in a generalized Hubbard model (with long-range 
interactions) on a lattice. The possible relation to a more realistic extended Hubbard Hamiltonian is discussed. 



PACS numbers; 71.10 Fd, 74.25.Ha, 74.25.Dw 

S.-C. Zhang recently conjectured that high- Tc cuprate 
compounds possess an approximate 5*0(5) symmetry. [Q. 
His theory aims to explain the proximity of supercon- 
ducting (SC) and antiferromagnetic (AF) phases in the 
phase diagram, and to account for the low-energy exci- 
tations as approximate 50(5) Goldstone modes. Anti- 
ferromagnetism and superconductivity are unified in one 
grand order parameter field (m^, m^, , Re^, Im^"), be- 
having as a 5-component vector, where the first three el- 
ements are Cartesian components of the staggered mag- 
netization and 5* is a spin-singlet SC order parameter 
(Here "Re^-" =^{^/ + '^^), etc.) 

In this picture, small symmetry-breaking terms tend to 
drive the system in a "superspin flop" between antifer- 
romagnctism and superconductivity, just as, in an mag- 
net with approximate 50(3) symmetry, competing spin- 
space anisotropics and external field can drive a "spin- 
flop" transition between magnetic states with order along 
the z axis and in the xy plane. 0] 

The 50(5) theory, while positing an intimate relation- 
ship between SC and AF order, does not imply that the 
pairing mechanism is AF fluctuations Rather, it 

quantifies the notion (also relevant to superfluid '^He) 
that there need not be a sharp difference between in- 
teractions mediated by magnetic and "charge" (number) 
fluctuations. I pass over Zhang's speciflc mechanism 
(whereby the system accommodates doping by switch- 
ing from the AF state to a symmetry-related SC state 
which has a different particle number), for 50(5) sym- 
metry can be valid even if another sort of perturbation 
is found responsible for the symmetry breaking and the 
AF-SC transition. 

The 41meV mode observed in spin-flip neutron scatter- 
ing on YBa2Cu307 |^ is interpreted as a Goldstone mode 
of 50(5) with a gap due to the symmetry-violating terms, 
analogous to the anisotropy gap in a spin-wave branch of 
the unixial magnet [0. These excitations are created by 
"tt'''" operators Q| (50(5) generators that mix magnetic 
and SC components), g They are charged bosons with 
the quantum numbers of "preformed" Cooper pairs, and 
presumably carry the current in the "normal" metal [Q ; 
it has been speculated that this explains the linear 
temperature dependence of the normal-state resistivity. 

To the extent that 50(5)-violating terms are small 
(as in Zhang's phase diagram "A" ||]), relations between 



AF and SC quantities are obviously predicted. For ex- 
ample, the Neel temperature T/v on one side ought to 
equal the SC Tc on the other side (the real ratio is 
5 : 1 in YBa2Cu307). Furthermore, when converted into 
the proper units, the tensors of superfluid density and 
AF spin stiffness should be equal; as should the order- 
parameter lengths (staggered moment and SC gap mag- 
nitude, respectively), and the interlayer couplings (inter- 
layer superexchange and intrinsic Josephson coupling, re- 
spectively.) An order-of-magnitude equality of the inter- 
layer couplings is indeed expected in the interlayer tun- 
neling picture Finally, the 50(5) Ginzburg-Landau 
theory predicts that vortices have magnetic cores [p],pT|; 
conversely, in analogy to the Bloch wall in the 50(3) 
magnet, it suggests that magnetic domain walls contain 
SC stripes, as proposed for other reasons by Emery and 
Kivelson 

Microscopic 50(5) symmetry — In this paper, using 
elementary notations, I implement a literal 50(5) sym- 
metry in a one-band lattice model, construct a Hamil- 
tonian with exact 50(5) symmetry, and finally consider 
whether a realistic Hamiltonian of an extended Hubbard 
form might approximate an 50(5) symmetric Hamil- 
tonian. Take a lattice with N sites (using periodic 
boundary conditions). Creation operators for the or- 
bitals on site x are cj.(x) for a = t, i- Let Q be the 
ordering wavevector of some two-sublattice AF state, so 
that e'*^ '' = ±1 at every site. The usual staggered- 
magnetization components are 



mM(x) ^ ie^Q-(c|(x)cT(x) -4(x)cx(x)), 



m^)(x)^e'Q-(4(x)ci(x)) 



(1) 



and m^''(x) = m^''(x)l'. The SC order parameter op- 
erator has the general form 'I'(x) = J^r r' "^(^^ -I- 
r)c|(x + r') (which allows different spatial symmetries 
depending on the form of the coefficients ip{r,r').) We 
seek a continuous, unitary operation that turns a com- 
ponent of m(x) into one of ^'(x), i.e. turns creation into 
annihilation operators: clearly it must be some form of 
Bogoliubov transformation. 

Indeed, a discrete 50(3) symmetry of this sort is al- 
ready known for the negative-?/ Hubbard model (is), for 
which the appropriate SC order parameter is \l/(x) — 
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C|(x)c|(x). One maps C|(x) e'^^'^Cl{^x.y (leav- 
ing C|(x) alone) which implies (vE'(x)'l', vl/(x), e**^'^n(x)) 
(to+(x), m_(x), m2(x)); here n(x) = c|(x)c|(x) + 

cj(x)c|(x). The only change induced in the Hubbard 
Hamiltonian is U —>■ —U; thus a hidden 5*0(3) symme- 
try relates SC order (^>) and charge-density-wave order 
{e^^'^n{x.)) in the limit of large negative U. 

To write the exact S0{5) symmetry transparently, and 
to ensure it in the order parameters and Hamiltonians, 
I use the duality |l6j between the "c" operators and an 
alternate set of canonically commuting operators, 



In real space, this says 

d,(x) = e-''«-'';^</p(r)c„(x + r) 



(2) 



(3) 



where ryk = X^r ^ (^(r). To make the symmetry ( ^ ) 
work, we will need the important conditions |15 



Vk+Q = -Vk 
?7-k = Vk 



(4a) 
(4b) 



for all k, which in real space say respectively that ifir) — 
for "even" r (meaning those connecting sites in the 
same sublattice) and that ip{r) = (/?(— r). Eq. ( |4a| ) implies 
{c^(x), dcr' (x')} = if X and x' are on the same sublattice 
(e.g. X = x'). Then the symmetry operation is just 

c;(x) = cos(0/2)c,(x) + sin(0/2)dL,(x) (5a) 
d'„ix) = - sin((/./2)cl,(x) + cos(0/2)d,(x) (5b) 

The symmetry ( pb| ) is generated by ^{n + n^), where 
TT = i'^-^[ci{x)di{K) — £ij^(x)c-f (x)]. To transform wave- 
functions, it is useful to know that the vacuum |0) tran- 



forms to J|j^(cos((/)/2) -I- ryt sin((/>/2)GLc 



,t 



-k+Q,i 



)|o). 



So that it will map exactly under (pa), the 5*0(5) stag- 
gered magnetization must be defined as 



m(x) 



1, 



m 



(x) 



m 



(d) 



(x) 



(6) 



where m'^''^ is (0) with "c"— Note this gives a sen- 
sible result for the Neel state: if m'^'^)(x) is up, then 
m(<i) (x) is down (since the d "orbital" on site x is a linear 
combination of "c" orbitals from the opposite sublattice). 
The SC order parameter is 



vI'(x)^e'Q--[ci(x)dT(x) 



■c^i(x)cT(x)] 



(7) 



Then 



mz'{x) = COS0 r7i^(x) 4- sin0 Re^(x) (8a) 
Re\l/'(x) = COS0 Re^I'(x) - sin0 r7i^(x) (8b) 



while the other three components are invariant. The 
50(5) rotation of the Neel state with (j) = tt/2 gives 



f7k4-fcLk|)|o) 



(9) 



This BCS state has no remnant of Fermi surface 
((c^Co-) = 1/2 throughout reciprocal space.) 

One could construct a total of six such rotations, each 
of which mixes one of the three components of m(x) 
with one of the two components of ^'(x). The other 
five could all be obtained by by combining (|5b| ) with the 
usual 50(3) rotations acting on the spin labels of c and d 
operators, plus the usual 5*0(2) = U{1) gauge symmetry 
changing their complex phases. (Zhang has discussed the 
algebra of 5*0(5) generators ||l|,|[.) 

For the square lattice, we must have Q — (tt, tt) so the 
hopping term dl^ ) will be 50(5)-invariant. This leaves 
much freedom to rj^, but the simplest choice is 



ryk = sign(cos — cos ky 



(10) 



This is inspired by the original and approximate 50(5) 
symmetry [|l| which had the same form but with coeffi- 
cients ?7k — » cos fca; — cosfcj,; recently Kohno j2j indepen- 
dently discovered the exact version (|l^). Comparison 
with (H) shows that (|l^) is essentially the Cooper pair 
wavefunction and has dx2-y2 pairing symmetry, consis- 
tent with strong experimental evidence in the cuprates 
|l7| . Interestingly, one other simple form would also sat- 
isfy the conditions (||): jyk = sign(cosfc2; -t- cosfcj^). That 
variant of 50(5), which entails "extended s-wave" pair- 
ing, appears free from internal contradictions (contrary 
to a suggestion in Ref. ||l|). 

The coefficients in (||) (Fourier transform of ([ic|)) are 
ip{x, y) = 4/[7r^(a;^ — j/^)] for x + y odd, zero for x + y 
even. For numerical and analytic explorations, it may 
also be helpful to have a one-dimensional toy realization 
of 50(5) symmetry. This is given by Q = tt and rjk = 
sign(cosfc), which gives (p{r) = 2(— l)(''"^-'-^/^/(7rr) for r 
odd, zero for r even. 

Microscopic Hamiltonian — Next I will produce an ar- 
tificial generalization of the Hubbard Hamiltonian which 
has exact 5*0(5) symmetry. The basic Hubbard model 
with particle/hole symmetry can be written 

nHubb^nhop + uJ2lH^)-'^)^ (11) 



n 



hop 



(-^)EE4(x)c.(x + u) =^ekcLck. (12) 



with u running over nearest neighbors, and Ck = 
{—t)(cos kx + cos ky). 

The minimal Hamiltonian including the terms in ( pT| ) 
is simply the 50(5) symmetrization of (pi]). The hopping 
term Tihop is already invariant imder all the 50(5) ro- 
tations (such as (8b)), provided that ek+Q = ^ Ck- That 
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is true in any bipartite lattice, if and only if Q describes 
the original Neel state with opposite spins orientations 
on nearest-neighbor sites. 

However, SO{5) symmetrization turns the number 
operator n(x) to something quite different, n*(x) = 
^ [n(x) — n^'*) (x)] . (I take the obvious definition n^'^^ (x) = 
(i|(x)(i|(x) + (i|(x)(i|(x).) In contrast to the usual num- 
ber operator, X^x^-^C^) = The n^'^\x) operator in- 
cludes terms |(^(r)pn(x + r), all on the opposite sublat- 
tice from x and largest for nearest neighbors, |r| = 1, 
as well as long range hopping between sites of the same 
sublattice. Thus the S'0(5)-symmetrized Hubbard model 
has a modified interaction term: 

=W/,op + C/'^^(n^(x))2 (13) 

X 

When we expand n'*(x)^, we get a variety of terms, which 
include interactions and hoppings (with diminishing co- 
efficients) to arbitrarily large distances. 

What is the ground state of (|l|)? If U" 0, at 
half-filling, it is the Fermi sea which manifestly possesses 
50(5) symmetry. The i — » limit of ( [T^ ) is more rel- 
evant and more challenging. Certainly a ground state 
is obtained at half filling by setting the even sublat- 
tice ferromagnetic in one direction and the odd sublat- 
tice ferromagnetic in any another direction, since n(x) — 
n^'^^(x) = 1 on every site. The two sublattice moments 
can be added to make a total angular momentum I tak- 
ing any value {0, 1, . . . , iV}. As was pointed out by 
Ref. Q|, by applying the tt and operators, as well 
as familiar spin-space rotations, each angular momen- 
tum is part of an 5*0(5) multiplet with a total degen- 
eracy + + 2) (2/ -I- 3)/6. This includes states 
with particle numbers differing from N by multiples of 
±2. The total degeneracy of this family of states is thus 
{N + l){N + 2)'^{N + 3)/l2, small compared to 2^ in the 
ordinary Hubbard model (|ll| ) with t — 0; however, con- 
ceivably this family does not exhaust the ground states. 

Now consider how a small t value splits these states 
in second-order perturbation theory. Among the states 
in which each sublattice is ferromagnetic aligned, this 
will have exactly the same effect as it does in the Hub- 
bard model. This suggests that the Neel state in fact 
approximates one of the ground states - and so must (||), 
the 50(5) rotation of the Neel state, since Ti. has 50(5) 
symmetry: thus I conjecture the t <^ V ground state 
has SO{5) broken symmetry. 

Group theory could be used to enumerate additional 
allowed terms in the Hamiltonian as in |2j; in particular, 
a bilinear coupling of the order parameter on neighbor- 
ing sites (50(5) symmetrization of the exchange inter- 
action). However, I have avoided this 50(5) t-J model 
analog. It could be derived (in the fashion I just out- 
lined, with I J| ~ t^/U') from the 50(5) Hubbard-model 
analog in the t W limit. But the most interesting 



phases of the standard t-J model occur in the regime of 
large J/t, which cannot be derived from any regime of 
the Hubbard model and the same thing may happen 
for the 50(5)-Hubbard model (|l|). 

Comparison to an extended Hubbard model — I now 
discuss how one might search for approximate 50(5) 
symmetry in some Hubbard-like model, such as 

He^t = HHubb + K,j, + + u) (14) 

xu 

In ( p^ ) Tijjop has the form of (p^), except that the coeffi- 
cient is t' and the displacement is u' running over second 
neighbors. The last term in ( p^ is a Coulomb repulsion 
between nearest neighbor sites. Comparison of measure- 
ments and calculations of electronic structure j2^ suggest 
that U/t > 4 and perhaps t' /t « —0.3 in cuprates. The 
aim is to find the point (s) in the parameter space of Tiext 
which make it closest to ([l^): can the parameters U, V, 
and t' of (|l|) be related to U" in (|l|)? 

We can very crudely guess at the U terms simply by 
retaining only the terms from ( [T^ ) of exactly this form. 
They come not only from n(x)^ but also from expanding 
of n('')(x)2. The resuh is U = iJ7"(l + IvWI*) = 
iC/«(l + l/9). 

We can estimate V in the same fashion, from the 
nearest-neighbor term in rt(x)7i*^'^'' (x) obtaining V ~ 
—W\(p{u.)\'^ where u is a nearest neighbor. Here the 
50(5) symmetry demands an attractive nearest-neighbor 
electron interaction, which is understandable: in the 
U'^ oo limit, only singly-occupied states can occur in 
a ground state, so only the AF states could be ground 
states. The SC state has a certain density of doubly- 
occupied and vacant sites, so an additional term is needed 
to equalize its energy with that of the AF state. Any 
pairing interaction might play the same role. 

Finally, no t' terms appear in (^^. In fact, any single- 
electron hopping term within the same sublattice vio- 
lates 50(5) symmetry and gets annihilated in the 50(5) 
symmetrization. Thus, although there are many quar- 
tic terms in (|lj) that do hop electrons between sites on 
the same sublattice, there is no 50(5) symmetric way of 
decoupling these terms to generate H'l^^p. (But in second- 
order perturbation theory, those quartic terms can gen- 
erate e.g. second-neighbor exchange interactions, just as 
the t' terms can.) 

Very recently, an extended Hubbard model has been 
diagonalized using an new interaction with double hop- 
pings, (l| Ex^W> where i^(x) = E.u4(x)c.(x + 
u) + C.C.. That is closer to dl^), since its terms have the 
same form as the largest terms (after those already men- 
tioned) in n(x)n''(x) and (n'^(x))^. This model has an 
apparently continuous AF/SC transition, so it may 
well realize 50(5) approximately. 

Of course, even at the 50(5) multicritical point in 
Zhang's picture, the microscopic Hamiltonian might have 
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no visible SO{5) symmetry; just as at the "spin-flop" 
point of an anisotropic magnet, a cancellation of terms 
favoring competing kinds of order might suffice, with the 
symmetry emerging only at long- wavelengths [Q. But 
if that length scale is much larger than the numerically 
tractable system size for Hubbard models, then direct 
numerical calculations on finite lattices (such as ) are 
too small to address the order parameter symmetry. 

Exact diagonalizations (e.g. [^) commonly study 
ground-state correlations, but their spatial decay is often 
inconclusive as a test of order due to small system size. 
Yet it is possible that the (excited) eigenstates show a 
well-defined structure characteristic of a particular sym- 
metry; this provided the convincing evidence for long- 
range order in the spin-1/2 triangular lattice AF pi| . 
(Very recently, Ref. [|l9| has pursued exactly such a pro- 
gram in exact diagonalizations of the t-J model). I 
suggest identifying an 50(5) multiplet numerically in 
a model with manifest SO{5) symmetry, |l9| ] and then 
following its evolution while the Hamiltonian is adiabat- 
ically modified to a more realistic model such as (|l4[). 

Conclusion — I have identified inklings of 5*0(5) sym- 
metry in popular existing models and exhibited the form 
an exact SO {5) symmetry could take in one- or two- di- 
mensional lattice models. The SO{5) symmetry in micro- 
scopic models is promising as a spur to the comparison or 
unification of competing models of high- Tc superconduc- 
tivity, and to improved understanding of extended Hub- 
bard models. However, I have not addressed the murkier 
issue of its application to the cuprates. 

Of the objections mounted so far to a possible SO{5) 
relationship between the actual AF and SC phases, one 
seems to be really inescapable: the Fermi surface If 
the SC metal shows a sharp drop in electron occupa- 
tion along a certain surface in reciprocal space, as found 
in angle-resolved photoemission experiments ||2^, then 
(see (H)) its AF image under SO{5) has a similar surface 
(shifted by Q). Apparently this AFM must be a spin- 
density- wave metal [Q. But the real AF phase of the 
cuprates is instead deemed to be a Mott insulator 
and its AF correlations are well modeled using a nearest 
neighbor exchange Hamiltonian [^ . 

I thank N.-G. Zhang for calling Ref. ^ to my attention, 
and S.-C. Zhang and S. Chakravarty for discussions. This 
work was supported by NSF grant DMR-9612304. 
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